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Abstract. We investigate advection dominated, trans- 
sonic accretion flows in the vicinity of a Kerr black hole. 
We take into account all relativistic effects in the dynam- 
ics of the flow and in the propagation of light. We as- 
sume the matter to be weakly magnetized and cool via 
the thermal synchrotron and Bremsstrahlung radiation. 
We include also the effects of Comptonization. We calcu- 
late the spectra of radiation as seen by observers located 
at different positions relative to the equatorial plane of 
the disk. The radiation emitted by the accreting matter is 
anisotropic and observers near the equatorial plane regis- 
ter a higher energy flux. This effect is more pronounced 
in the case of slowly rotating black holes. We calculate 
also the shape of a hypothetical gamma line, which may 
be produced by the thermonuclear reactions in the inner 
part of the flow. The line is strongly broadened, but the 
fact that the flow is quasi-spherical removes the two-peak 
shape of the line seen in the spectra emitted from thin, 
Keplerian disks. The kinematics of the advection domi- 
nated flows is not unique (as opposed to Keplerian disks 
or spherical free-fall) and it would probably be difficult 
to find strong limits on source models using the spectral 
observations. 
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1. Introduction 

The advection dominated accretion flows (ADAFs) have 
recently become an interesting alternative to the tradi- 
tional thin accretion disks (Shakura & Sunyaev 1973; 
Novikov & Thorne 1973). The disk model with a two 
temperature, optically thin and hot inner region has been 
proposed by Shapiro, Lightman & Eardley (1976). The 
spherical accretion flow, which cannot cool, because the 
radiation is trapped and advected with matter, has been 



considered by Katz (1977). Spherical accretion with ad- 
vection has been investigated by Begelman (1978, 1979). 
The role of advection in disk-like accretion has been stud- 
ied by Begelman & Meier (1982), Abramowicz, Lasota & 
Xu (1986) and Abramowicz et al. (1988) for configurations 
with substantial optical and geometrical thickness. These 
models are sometimes called slim accretion disks. The ad- 
vection of heat makes the energy conservation nonlocal, as 
opposed to thin disks, and thus influences their structure, 
but due to the substantial optical thickness, has no promi- 
nent influence on their spectrum. Another branch of ac- 
cretion flows has been discovered by Narayan & Yi (1995) 
and Abramowicz et al. (1995). Flows belonging to this 
class are extremely optically thin and have low accretion 
rate. Because of the low radiation efficiency, the gas can 
heat up to near virial temperatures. For this reason the 
flow becomes quasi-spherical and pressure gradients play 
an important role in the gas dynamics. The radial and az- 
imuthal velocities can have the same order of magnitude. 
These are the main characteristics of ADAFs. Because of 
their high temperatures ADAFs produce X-rays in a nat- 
ural way. On the other hand the effectiveness of accretion 
is low, because almost all the thermal energy goes under 
the horizon with accreting matter. These properties make 
ADAFs possible explanation for some X-ray sources and 
such models have been proposed (Narayan 1996; Lasota et 
al. 1996; Narayan, McClintock & Yi 1996; see also Rees, 
1982). 

The general properties of the advective flows have 
been obtained with the help of Newtonian or pseudo- 
newtonian (Paczyhski & Wiita 1980) description of the 
background gravitational field. Lasota (1994), Abramow- 
icz et al. (1996, hereafter ACGL) and Abramowicz, Lanza 
& Percival (1997, hereafter ALP) have obtained a self- 
consistent system of equations describing the accretion 
flows in the Kerr spacetime. All these papers use the verti- 
cal averaging to obtain a set of ordinary differential equa- 
tions to describe the flow. The form of equations resembles 
the equations for an accretion disk, but some properties 
of the quasi-sphericity are taken into account (see ALP). 
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The properties of ADAFs in fully relativistic treatment 
have been also investigated by Pcitz & Appl (1996). All 
these papers investigate mainly the flow dynamics and 
are concentrated on its properties and their dependence 
on parameters such as viscosity parameter a and the di- 
mcnsionless accretion rate rh. 

We are interested in the radiation processes and spec- 
tra emitted by ADAFs in the Kerr spacetime. The main 
reason of our investigation is to get a fully self-consistent, 
relativistic description of the radiation from the advec- 
tive flows. The gas has the highest temperatures near the 
horizon, so the combined effects of ray deflection, Dopplcr 
shifts and gravitational redshift may influence the high 
energy end of the spectrum, at least quantitatively. Since 
the flow is optically thin, the observer may have a chance 
to see the matter on the other side of the black hole ap- 
proaching the telescope with relativistic speed; in some 
cases this should lead to the beaming of radiation and sub- 
stantial increase of the surface brightness. The possibility 
and importance of such effects should be investigated. 

Since we are mostly interested in the propagation of 
light, we are not going to investigate the hydrodynamics 
of the flow in detail. To get rough solutions of the hydro- 
dynamic equations we use an approximate treatment, ne- 
glecting the cooling. Thus the most time consuming part 
of the calculations is decoupled from the equations of mo- 
tion. For a simple equation of state the structure of the 
flow (i.e. the profile of the angular momentum, radial ve- 
locity of matter and the speed of sound) depends only on 
the viscosity parameter a and the values of the specific 
angular momentum and energy of matter at the horizon. 
The matter density depends on the above variables and 
scales with the accretion rate. Using the same solution of 
the equations of motion one is able to get models with 
different accretion rate. Inspection of the solutions of the 
hydrodynamic equations of Chen, Abramowicz & Lasota 
(1997), which include cooling, shows that dependence on 
the accretion rate is weak, at least in the inner part of the 
flows. 

In the next Section we present the system of equations, 
which we use to model the accretion flow and rays trac- 
ing. We also describe the approach to radiation processes 
within the flow. In Sec. 3 we present the results of calcula- 
tions, showing the brightness profiles of our sources, their 
spectra, and luminosity dependence on the observer incli- 
nation angle. The discussion follows in the last Section. 

2. The model 

2.1. The accretion flow 

We investigate the stationary flow of matter and the prop- 
agation of light in the gravitational field of a rotating Kerr 
black hole using the Boyer-Lindquist coordinates t, <f), r, 9 
and the metric components g a b(r,9) as given by Bardeen 
(1973). We follow the (—,+,+,+) signature convention. 



When dealing with the flow we use the metric tensor with 
components evaluated at the equatorial plane (9 = tt/2). 
We use geometrical units, so the speed of light c = 1 
and the mass of the hole M = 1. (Later we use c for the 
speed of sound.) The Kerr parameter a (0 < a < 1) gives 
the black hole angular momentum in geometrical units. 
We use the Einstein summation convention, where needed 
and a semicolon for the covariant derivative. The normal- 
ization of four velocity with the chosen metric signature 
reads u a u a = — 1. 

The system of equations we use follows in general 
ACGL. Since we are neglecting the cooling processes at 
this stage, and treat the accretion flow in the disk ap- 
proximation, two components of velocity and the speed 
of sound c given as functions of radius, fully describe the 
dynamics. We use the angular velocity SI and the physi- 
cal radial velocity V as measured by locally nonrotating 
observers as main kinematic variables (compare ACGL). 
After the vertical averaging the velocity perpendicular to 
the equatorial plane is neglected (u e = 0) and other com- 
ponents of the four velocity are given as: 



u* = 



1 



y/T=V*y/-gtt - 2f! 3t0 - W 



flu 1 



V 



fg~^ Vl - V 2 



(1) 



(2) 



(3) 



We shall later use the two temperature plasma to de- 
scribe the cooling processes. In the hottest and most inter- 
esting part of the flow the pressure of the nonrelativistic 
ions, p — nikTi, where rii is ion concentration, Ti their 
temperature and k the Boltzmann constant, dominates. 
The isothermal sound speed can be defined as 



m 

m 



(4) 



where m is the average ion mass, so one has p — poc 2 , 
where p is the rest mass density. The specific enthalpy p 
is given as: 

H = = 1 + -& 

Po 2 



(5) 



where e is the total (rest mass plus thermal) energy den- 
sity. The above formula implicitly assumes that the matter 
can be treated as the ideal, nonrelativistic gas. 

The mass conservation equation relates the mass ac- 
cretion rate M, with the radial velocity and the surface 
rest mass density £ (ACGL): 



M = -2ttSv / A 



V 



(6) 
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where A(r) = r 2 — 2r + a 2 is one of the functions defining 
the metric. On the horizon one has A(r^) = and V = — 1 
(see ACGL), but the ratio A/(l — V 2 ) remains finite and 
is related to the derivative of radial velocity at r = r^. 
Thus the surface mass density remains finite everywhere 
and using mass conservation one can substitute for £ in 
other equations. 

The rest mass density of matter can also be replaced 
in the equations with the help of the so called vertical 
hydrostatic equilibrium equation, which reads: 



space time. The projection tensor = 5^ + u a Ub and its 
3 d trace = 2. This is why we use the factor 1/2 instead 
of conventional 1 /3 in front of the term including the ve- 
locity expansion 9 = u^ a in the definition of the shear 
tensor. 

The projected equation of motion, h h r T£. a = takes 
the form: 



V 2 



c 2 /pV c 2 c' _ 1 c 2 A' 



1 - V 2 V +2 p c 2 fj, A 



dp 



(7) 



where g jZ represents the gravitational tidal forces as mea- 
sured at the equatorial plane by an observer comoving 
with the fluid, and z is his coordinate in this direction. 
We use the prescription of ALP for the tidal forces, so the 
expression remains finite at the horizon: 



9.z 



-a 2 {u 2 t ~l) 



£ 1 

r-4 



(8) 



For an optically thin medium we assume the temperature 
to be approximately constant in the vertical direction. The 
integration of of the equation gives the density distribu- 
tion: 



1 z 2 

Po (r, z) = poM)cxp ( ~2^2 



(9) 



+-(i-v 2 )(u t ) 2 (g' tt + 2ng>, + n 



2, J 



(13) 



where the primes denote the derivatives relative to r. We 
neglect the terms related to the viscosity in this equation, 
which is the usual simplification (ACGL, Peitz & Appl 
1997). 

Projecting the equation of motion into the comoving 
reference frame, u b T£. a = leads to the energy conser- 
vation equation in the form (Novikov & Thornc 1973, 
ACGL): 



p Tu a V a S = 2na 2 



(14) 



where £ is the entropy of the fluid per unit mass and 
the derivative is in fact relative to the proper time. In 
the stationary, axially symmetric case u a V a — u r d r . The 
entropy of an ideal gas is given as: 



where 
H 2 = 



2 4 



(10) 



S = — (| In T - In p ) + const (15) 

771 A 

For the viscosity in the a-disk approximation one has: 



is the scale height. The form of the above expression 
shows, that c 2 / ' p has the physical meaning of the square of 
the speed of sound. Integration of the density distribution 
over the disk thickness gives £ = \Z~2Trp H. Knowing the 
surface mass density and the scale height one can substi- 
tute for the rest mass density. 

We are neglecting transport of heat within the fluid, 
but we do include the processes due to the shear viscosity, 
so the energy momentum tensor has the form (Novikov & 
Thorne 1973): 



T b a = (e+p)u a u b +p6Z-2no-Z 



(11) 



where 5% is the Kronecker delta, n the dynamical viscosity, 



and al the shear tensor: 

O-ab = ^{Ua;iK + U b:i h l a ) - ^Qh ab 



(12) 



The equations of motion for the disk are written under the 
assumption that the velocity in the vertical direction can 
be neglected and all variables are evaluated in the equato- 
rial plane. The motion is effectively two dimensional and 
all tensors should be treated as acting in three dimensional 



V = ap —H = apo — — 



(16) 



where we use the expression for the sound speed with the 
proper physical meaning. The equation for the entropy 
gradient can now be expressed as: 



V c! H' 1A' 2ar i /(l - V 2 ) ., 

T~ ' ''7 77 ° 



AV 2 



(17) 



where the radial derivative of the scale height is still to 
be expressed via the gradients of the angular and radial 
velocities and the speed of sound. 

The symmetries of the Kerr metric imply the energy 
and angular momentum conservation equations for the 
fluid in the form T?. a = 0, T?. n = 0. After integrating 
over r and vertically averaging one has: 



M c 2 r 2 
_ (w -M°<) + 2arE — 







- ( ^-A?) + 2^— 



(18) 



(19) 
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where the superscript denotes variables as measured on 
the horizon. We define the specific energy of the fluid as 
£ = —jj,Ut and its specific angular momentum as L — \xu$. 
We also use the kinematic specific angular momentum I = 
—u^/ut- Combining the upper equation multiplied by £ 
with the lower equation multiplied by L we get: 



ac 2 r 3 



AV 2 



1-V 2 



££o(t-£o) 



(20) 



The combination of the shear tensor components, which 
we use, contains the angular velocity gradient but not the 
radial velocity gradient: 



Mv*r v 2 [£g' tt + (i + ne)g^ + n g ' ct)i 
1 I '- " 



(1 - V 2 ) 2 (l - Q.I) 



(21) 



Thus the equation for O is obtained. It has a boundary 
condition at the horizon = w, where u> = —g^jg^ is 
the angular velocity of the dragging of inertial frames. 

We solve the equations for V, c and by iterations. 
We start from a solution with subkeplerian rotation, which 
asymptotically (far from the hole) becomes a selfsimilar 
solution as described by Narayan & Yi (1994). We postu- 
late the position of the sonic radius r s and the value of the 
angular momentum at the horizon £q and, using the angu- 
lar velocity distribution from the initial guess (and later 
on from the previous iteration), we solve the equations for 
the radial velocity and the speed of sound. These equations 
are easy to integrate if one starts from the sonic point. For 
a given angular velocity distribution, in particular know- 
ing its value at the sonic point Sl we find the values of 
the other two variables V s , c s from the regularity condi- 
tions. The boundary condition on the horizon (V = — 1) 
is always met. The value of specific energy on the horizon 
results from the solution and is iterated. Next we solve 
the equation for the angular velocity starting from the 
horizon outward, which gives a new angular momentum 
distribution. The expression for the square of the shear 
{a 2 ) contains derivatives of the angular and radial veloc- 
ities, and is rather complicated, including the derivatives 
in the second power. Using Eq.(12) we calculate the shear 
tensor from the "previous iteration" and substitute it into 
Eq.(17). 

Not all postulated (r Sl £ ) pairs lead to solutions which 
can be continued from the horizon to infinity, but the it- 
erations just described are so fast that checking many of 
them is easy. For a given angular momentum of the black 
hole we find the best parameters, for which the solution 
can be smoothly joined with the selfsimilar solution far 
away. 

For the purpose of investigating the emitted spectra we 
choose examples of ADAF solutions in Kerr metric with 



> - 



3 
G 




Fig. 1. Sample models of ADAFs. On the upper plot we show 
the radial velocity of the flow V (negative) and the sound speed 
c* = c/ \/Jl as functions of 1/r for our three models in the Kerr 
metric with angular momentum a — (dotted line), a = 0.5 
(dashed line) and a — 0.9 (solid line). The locations of the 
sonic points are shown as large dots on the lines. On the lower 
plot we show the angular velocity distribution f2 — lo, using the 
same conventions. 



angular momentum a = 0, 0.5 and 0.9. We use the vis- 
cosity parameter a = 0.1. For the Schwarzschild case the 
position of sonic radius and the value of the angular mo- 
mentum on the horizon are in good agreement with Chen 
et al. (1997), who use the pseudo-Newtonian potential to 
describe gravity. The profiles of the sound speed and ra- 
dial and angular velocities for our solutions are shown on 
Fig.l. 

2.2. Radiation processes in the disk 

We use two temperature plasma to describe radiation pro- 
cesses and local spectra in the disk. For given ion tem- 
perature Ti, scale height H and surface mass density S 
(Section 2.1.) we determine electron temperature by tak- 
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ing into account the detailed balance of their heating and 
cooling. We divide the disk into the coaxial cylindrical 
layers and treat them separately. We imagine each of the 
annuli to be cut from a plane parallel layer with con- 
stant temperature and the density concentrated toward 
the midplane according to Eq(9). We follow here the de- 
scription of heating and cooling processes presented by 
Narayan & Yi (1995) and also by Esin ct al. (1996) [and 
references therein] . They consider the energy transfer from 
ions to electrons via Coulomb collisions and cooling of elec- 
trons by synchrotron radiation, inverse Compton process 
and brcmsstrahlung emission. They ignore the presence of 
e+e~ pairs but as Bjornsson et al. (1996) showed, the role 
of e + e~ pairs in the structure of ADAF is not significant. 

For calculated electron temperature we can model the 
local continuum emission per unit area from surface of 
each annulus. We assume that molecular weights of ions 
and electrons are fii — 1.29 and p e = 1.18, and the 
density of magnetic field is parameterized by the ratio 
of gas pressure to total pressure f3. In our calculations 
we ignore the radiation pressure, so the total pressure 
is the sum of gas pressure and magnetic field pressure. 
We use the approximate formulae of Esin et al. (1996) 
to calculate synchrotron spectrum and formulae of Jones 
(1968) to model the inverse Compton scattering of syn- 
chrotron photons. The method of approximate compu- 
tation of brcmsstrahlung spectrum is given in Svensson 
(1982) [and references therein]. We include the frequency 
averaged comptonization of brcmsstrahlung photons in 
the energy balance, but we neglect this process in the cal- 
culation of the spectrum, since it is very time consuming 
and its contribution to the resultant spectrum is not sig- 
nificant. 

Figure 2 shows the continuum emission calculated 
for four annuli of radii r/r g = 1.5; 10; 100; 1000, (where 
r g = GM/c 2 is the gravitational radius) and for the fol- 
lowing parameters: M/M Q = 3.6 • 10 7 , M/M Edd = 0.016, 
corresponding to 1.5 • 10 _3 M Q yr _1 , a = 0.1 and (3 = 0.95. 
Such parameters have been chosen by Lasota et al. (1996) 
to model the spectral features of the nucleus of active 
galaxy NGC 4258, which they treat as a superposition of 
ADAF and a standard accretion disk. We show the hottest 
of our models, which has a — 0.9. 

The importance of the Compton scattering depends on 
the total optical depth of the disk and, strictly speaking, 
is not a local process. In the following part we neglect this 
fact. The bremsstrahlung emission is a two body process 
and, under the assumption of constant temperature in ver- 
tical direction, we assume its rate to be proportional to the 
square of the density (~ p 2 ). Similarly the synchrotron 
emission, being proportional to the product of electron 
density and magnetic pressure, which has a constant ra- 
tio to the gas pressure and density, is also proportional to 
the square of the electron concentration. The rate of the 
Compton process depends on the electron concentration 
and the density of the radiation field. The latter one de- 




lg(") 

Fig. 2. Continuum emission per unit area from four annuli of 
radii r/r g — 1.5; 10; 100; 1000 (reading from highest line down) 
for parameters: M/M Q = 3.6 • 10 7 , M/M Edd = 0.016, corre- 
sponding to 1.5 • 10~ 3 M o yr~\ a = 0.1, (3 = 0.95, and a = 0.9. 
The left peak of each curve is due to synchrotron radiation, the 
right one is due to bremsstrahlung emission and two peaks be- 
tween them are due to the Compton scattering of synchrotron 
photons. 

pends on the optical depth for scattering measured to the 
disk surface. Thus for this process the proportionality to 
the square of the matter density does not hold, but the 
decline toward the surface is also present and should have 
similar character. We use the following, simplified formula 
for the volume emissivity from the disk: 

[2 F(v) ( z 2 \ 

The flux here is the total energy emitted from the unit 
surface of the disk as usually used in this context. The 
emissivity j(y) is defined per one steradian and hence the 
integration of the above formula from the equatorial plane 
to infinity would give F(u)/n. 

Our formalism described in the next subsection allows 
for the calculation of the source spectrum regardless of 
the nature of the radiation processes in the disk as long 
as the matter remains transparent to the radiation. We 
consider also the shape of the emission line, which may be 
used for testing the velocity field in the flow. We follow the 
idea of Narayan & Yi (1996), who consider the spallation 
process leading to the lithium and beryllium production. 
As shown by Kozlovsky & Ramaty (1974) some of the nu- 
clei are produced in the excited state and emit 7-ray lines 
at 478 and 431 keV. The emissivity is proportional to the 
rate of production. Following closely Narayan & Yi (1996) 
we assume that there is a threshold energy 8.5 MeV for 
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the a-a reaction to occur and we assume the crossection 
to be inversely proportional to the energy of interacting 
particles a <~ E~ 2 <~ c~ 4 (for c > c t hroshoid)- We use here 
the typical thermal energy of the a particles instead of 
energy distribution averaging. The relative velocity of in- 
teracting particles is proportional to the thermal energy, 
and their concentration to the matter density. The width 
of the line results from the thermal motion of 7 Li ions; 
they are slower than other particles because of the higher 
mass, and including the \J~2 factor in the definition of the 
line width, we get: 



Ais 



1.1c 



(23) 



where vq is the line frequency and the mean ion mass 
1.29™^ of the matter containing 70% of hydrogen has 
been assumed in the calculation. Taking all the factors 
into account we obtain the formula for the line volume 
emissivity up to a constant factor: 



2 CXp 



\ (A^) 2 ) 



(24) 



where the density is obtained from Eq.(9), and the formula 
is valid only for the sound velocity over the threshold. 

2.3. Geometrical optics 

We use the backward ray shooting method (Luminet 1979; 
Jaroszyriski 1993) to obtain the spectrum of the electro- 
magnetic radiation from the configuration. The intensity 
of radiation I(fo) coming to a far observer at the frequency 
vq along any ray is given as: 



J ray ^ 



(25) 



where is the emissivity per unit volume, and dl pr0 p 
is the proper length differential as measured in the frame 
comoving with matter. The emission frequency v is given 
as: 



U Pa 

Po 



(26) 



where u a is the matter four velocity, p a is the photon 
four momentum, and po - the photon energy in BL frame, 
which is conserved. This is the usual redshift (or blueshift) 
formula. 

We find the null geodesies in the Kerr metric using the 
quadratures of Hamilton- Jacobi equation (Bardeen 1973). 
This method is faster, and more importantly, has a bet- 
ter accuracy as compared with the solution of ordinary 
differential equations along the trajectories. 

Following many rays coming to an observer from the 
vicinity of the hole makes it possible to produce a map of 
the surface brightness of the source. Since the hottest part 
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Fig. 3. Bolometric brightness profiles as seen by observers at 
the angles 6 = 5° to 85°, bottom to top. The radius is in 
geometrical units, the brightness in c.g.s. units. The three plots 
correspond to the three models as indicated. 



of the source is close to the hole, we probe the radiation 
intensity in rings centered on the hole with logarithmic 
distribution of radii. Thus the central parts are investi- 
gated with greater accuracy. A useful characteristic of the 
source is its radial surface brightness profile, as seen by a 
distant observer, obtained by integration over frequency 
and azimuthally averaging: 



i 



(27) 
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where £ and are the observer cylindrical coordinates on 
the sky. The length corresponding to the angle £, measured 
at the source, is R = £D, where D is the distance to the 
source. The flux of radiation from the source as measured 
by an observer is given as: 



and implicitly depends on the observer location. 



(28) 



3. Results 

We have made ray shooting calculations for observers at 
inclination angles 9 = 5°, 30°, 40°, 50°, 60°, 65°, 70°, 75°, 
80° and 85° as measured relative to the rotation axis, for 
the three models of ADAF in three background metrics 
with the Kerr parameter values a — 0, 0.5 and 0.9. In 
ray tracing we include only the radiation from a sphere of 
the radius r = 10 3 r 9 . Wc do not include any other source 
components, as a thin Keplerian disk or a dust torus. 

In Fig. 3 we show the source surface brightness profiles 
as seen by observers at different inclinations. The rela- 
tively low surface brightness in the centre corresponds to 
rays, which traverse only the matter in front of the hori- 
zon. The densest and hottest part of the matter (close 
to the horizon) is receding from the observer, which also 
works as to decrease the observed brightness. The abrupt 
jump to higher surface brightness corresponds to rays 
which pass close to the hole. These rays spend relatively 
long time in the innermost part of the flow due to the 
bending. Some of them may be tangent to matter trajec- 
tories in isolated points. Dopplcr boost and relative in- 
crease in optical depth (which remains much lower than 
unity) both work as to increase the surface brightness. The 
dependence on the inclination angle is better pronounced 
for the low angular momentum black holes. The surface 
brightness is higher, when the source is watched from the 
vicinity of the equatorial plane. 

The bolometric flux is also higher for observers close 
to the equatorial plane. As can be seen in Fig. 4, the effect 
is moderate for the fast rotating hole (about 70% increase 
in flux measured at the equator as compared to the pole). 
For slowly rotating holes (a = 0.5 or a = 0) the source 
seen from the equatorial plane seems to be brighter by a 
factor 4. 

The shape of the continuum spectra from our models is 
not considerably changed by the relativistic effects in the 
ray propagation. As seen on Fig. 5 the spectra are shifted 
vertically, when the observer approaches the equatorial 
plane. The simultaneous shift to higher frequencies is very 
weak. 

The line profiles are obtained with the help of the 
same calculation scheme, but the volume emissivity of 
matter given in Eq.(22) is replaced by Eq.(24). The lines 
are appreciably broadened. The mechanism of line forma- 
tion which we use (spallation reactions leading to the pro- 
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Fig. 4. The dependence of the observed total bolometric flux 
from the source on the observer inclination relative to the axis 
of rotation. We use the bolometric flux from our Schwarzschild 
model (a — 0) as seen by an observer on the rotation axis as 
the unit of flux. The dotted line corresponds to the nonrotating 
hole (a = 0), the dashed line to the case a = 0.5 and the solid 
line to a — 0.9. 

duction of lithium and beryllium nuclei in excited state), 
gives the lines which are broadened by thermal motion 
of the hot gas. The effects in light propagation (Doppler 
and gravitational shifts) are superimposed on the natu- 
ral line broadening. As can be seen on the plots, these 
combined effects produce profiles without the character- 
istic two-peak structure as in the case of Keplerian disk 
(Hameury, Marck & Pelat 1994). Such two-peak structure 
is present in the spectra of the maximal surface bright- 
ness rings (compare Fig. 3), but the light from these parts 
of the sources cannot be observed separately by a distant 
observer. Superposition of the spectra from all the rings 
removes the effect, but some traces of it can be seen in the 
spectra of our a = and to lesser extent a = 0.5 models. 

4. Discussion and conclusions 

Our system of equations describing the relativistic advec- 
tion dominated accretion flow onto a Kerr black hole is 
based on the ACGL and ALP papers but there are few 
differences. We use all the terms related to the radial and 
azimuthal motion of matter in the expression for shear, 
which is responsible for heat dissipation. We also combine 
the equations for energy and angular momentum conser- 
vation to get the equation for the angular momentum gra- 
dient, which does not contain gradients of radial velocity. 
(ACGL neglect such terms). We neglect the viscous forces 
in the equation for radial momentum balance; this ap- 
proximation is also used by ACGL). We take into account 
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Fig. 5. The spectra of our models for observers at inclina- Fig. 6. The hypothetical line profiles emitted by our models for 
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(thick lines). The three panels correspond to the three models lines), and 85° (thick lines). The three panels correspond to the 
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the contribution of the thermal energy to the total energy- 
density introducing the specific enthalpy fi into our equa- 
tions. Our description is however limited to the equation 
of state for an ideal, nonrelativistic gas. This approxima- 
tion is reasonable, since the thermal energies of ions and 
electrons do not exceed their rest energies. The contribu- 
tion to pressure of the magnetic field is low (5%) and does 
not limit the validity of the ideal gas approximation. 

The most serious problem in our approach (and in 
other descriptions of ADAF in the frame of vertically aver- 
aged disk structure) is the fact, that the flow is not limited 



to the close vicinity of the equatorial plane. We obtain the 
three dimensional description of the flow assuming that 
the flow is quasi-spherical (we neglect the component u e 
of the velocity) and that the density falls in the vertical di- 
rection as in isothermal atmosphere, with a constant scale 
height depending only on radius. 

The influence of the relativistic effects in light propa- 
gation on the shape of the continuum spectra of advection 
dominated accretion flows seems to be marginal. This is 
due to the fact that the emitted spectra are very broad 
and the small frequency shifts arc hard to distinguish. 
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The most important effect we hnd is the anisotropy 
of the radiation emitted from the disk, which favors ob- 
servers near the equatorial plane. Observers close to this 
plane can see part of the matter as approaching them with 
the relativistic speed due to the combined effects of rota- 
tion and fast radial fall. This effect can be stronger that 
the occurrence of hot spots on the surface of thin Keple- 
rian disk in the places approaching the observer with the 
highest velocity (Luminet 1979; Jaroszyhski, Wambsganss 
& Paczyhski 1992). Since matter is transparent one can 
see the beaming of radiation regardless of the place of its 
emission. The fact, that the effect is stronger in our models 
with slowly rotating holes is due to the different angular 
momentum distributions in the accreting matter. Also the 
fact, that the flow on the fast rotating hole (a = 0.9) has 
a higher ion temperature (see the plots of the sound speed 
on Fig.l) has an impact. Since the electron temperature 
does not grow with ion temperature indefinitely but sat- 
urates at about 10 10 K, it makes the central, fast rotating 
part of the flow relatively less important as a source of ra- 
diation and relativistic effects become weaker than usually 
expected for a fast rotating hole. 

The shape of the hypothetical 7-ray line, which we 
calculate, does not have a two-peak structure due to the 
superposition of photons coming from the different parts 
of the flow. The specific production mechanism of the line, 
which we use, also acts as to make the relativistic effects 
less pronounced. The crossection for the lithium and beryl- 
lium production decreases with the ion temperature. The 
hotter flow on the Kerr hole with a — 0.9 produces a 
stronger line, since the volume with the temperature above 
the threshold is larger, but in the same time, the central 
region is much less important. 

In the conclusion we state, that the kinematics of ad- 
vectively dominated accretion flows is difficult to investi- 
gate on the basis of spectral observations, since the pos- 
sible pronounced spectral features are lacking. Also the 
parameter space of the possible ADAFs is much larger 
than in the case of Keplerian disks, which makes the task 
even harder. 
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